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Modified gravity models are subject to a number of consistency requirements wliicli restrict the 
form that the function F{R) can take. We study a particular class of F{R) functions which satisfy 
various constraints that have been found in the literature. These models have a late time accelerating 
epoch, and an acceptable matter era. We calculate the Friedmann equation for our models, and 
show that in order to satisfy the constraints we impose, they must mimic General Relativity plus A 
throughout the cosmic history, with exponentially suppressed corrections. We also find that the free 
parameters in our model must be fine tuned to obtain an acceptable late time accelerating phase. 
We discuss the generality of this conclusion. 

PACS numbers: 



In recent years, experimental cosmology has provided strong evidence that the Universe is currently undergoing a 
phase of accelerated expansion. The physical phenomena responsible for this expansion has proved elusive, and there 
are numerous theories that attempt to explain the current epoch of the Universe. Typically, the observed acceleration 
is attributed to a new energy component, dark energy, which dominates at late times. Numerous models exist of this 
form, the simplest of which include introducing a cosmological constant A, or alternatively a scalar field slowly rolling 
down a potential. 

Although dark energy is undoubtedly the most popular explanation of the current epoch of the Universe, it is 
not the only way to obtain late time acceleration. In this paper we consider a class of models where gravity is 
modified at large scales, in such a way that the late-time expansion of the Universe could arise naturally from the 
new gravitational field equations. These modified gravity models have been considered by numerous authors, see for 



. example [l|, 1, S 1 i, i, & i i, 0, [n|, [ll, 0, Q & E [ll 



Modifying gravity in a manner consistent with experimental data has proved difficult. General Relativity (GR) is a 
very robust and well tested theory, and it has been found that even slight modifications often lead to instabilities, such 
as the propagation of ghosts [ly, or in otherwise stable matter sources 0,[2']. When introducing additional terms 
into the gravitational action, we must be careful to respect the success of GR in both the low and high curvature 
regimes, to ensure that any potentially new model agrees with known observational tests of gravity. 

In the literature, there are numerous simple conditions that exist which restrict the form that a modified gravity 
model may take. In this paper, we combine all of these conditions, and in doing so find that the modified gravity 
function F{R) is well constrained, even before considering stringent experimental tests. The philosophy adopted 
throughout is that before attempting to reconcile any potential modified gravity model with experimental data, that 
modified gravity model must first satisfy these consistency requirements. 

We will take a trial F{R) model which satisfies all of the existing constraints. We find that our model exhibits 
a late-time accelerating solution, in the absence of a cosmological constant. We then consider the cosmology of 
our model, and obtain a modified Friedmann equation. We show that to satisfy all constraints, and to obtain an 
acceptable Newtonian limit for small R, our model must act like GR with highly suppressed corrections throughout 
its cosmological history. In addition, wc find that although our model has no true cosmological constant, we must 
still fine tune a particular combination of the free parameters in F{R) to obtain an acceptably small curvature in the 
current accelerating epoch. We conjecture that our results are generic features of models which satisfy the consistency 
requirements. 



Before we consider a specific F{R) model, we first briefly review the fleld equations in modified gravity theories. 
Throughout this paper, we will study the action 



I. INTRODUCTION 



II. RELEVANT FIELD EQUATIONS IN THE JORDAN AND EINSTEIN FRAMES 
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where F{R) is some function of the Ricci scalar R, M is the mass scale of our model, and S'm the matter action. The 
field equations obtained by varying the action ([T]) with respect to the metric are 



F'iR)R^,, - \F{R)g^, - V^V^F'iR) + g^.UF'iR) = ^T^,, 



(2) 



where T^^ is the energy momentum tensor of any matter present. Taking the trace of ^ gives the more useful 
equation 



Q{R) + 2.UF'{R) = ^, 



(3) 



where the function Q{R) is given by Q{R) = RF'(R) — 2F{R). Vacuum solutions to the field equations satisfy 
R = Ro = const when T = 0, and hence Q{Ro) = 0. 

In addition to the above field equations in the Jordan frame, we will also want to consider the potential of the 
scalar field in the Einstein frame. To transform to the Einstein frame, we first write ^ as 



S 



^gd^x[F'{c^) {R-ct>) + Fm 



(4) 



where £„! is the matter Lagrange density. The actions (jl]) and ([T]) are equivalent if F"{R) ^ 0. We then use the 
conformal transform g^i, = F'((j))gfj_,y, which gives 



(5) 



R = F'{cf,)[R + 3a\nF' -^^Va.F'S/"F' 



where R and R are the Ricci scalars associated with the metrics g^j^ and respectively. Using this in ^ gives the 
following Einstein frame action 



= y \J —gd'^x 
where the Einstein frame potential is given by 



(6) 



^^^^ ^ — \ F'{^Y 
Finally, by making the field redefinition F'[(^) — exp(crV2/\/3AP), we obtain 



(7) 



gd'^x 



AP ~ 1 / - \2 



1 



with potential 



(8) 



„cr\/2/V3M2 



-2(T\/2/\/3M2 



(9) 



In many instances, it is preferable to study modified gravity in the Einstein frame, since the field equations are much 
simpler. Additionally, the potential ^ is useful since its critical points correspond to the vacuum solutions Q{Rq) = 
in the Jordan frame. We stress that the two frames are globally equivalent only if F"{R) ^ 0. 
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III. REVIEW OF CONSTRAINTS ON F{R) MODELS 

Here we briefly review the constraints in the literature imposed on the form of the function F{R). We begin by 
demanding F'{R) > and F" ^ for all R. The first of these is to ensure stability; it has been shown [l^ that if 
F'{R) < for any R then the model will possess ghost instabilities in this region. The second condition F"{R) ^ 
is imposed since it has been found that models which have F"{R) < will generically possess instabilities, when we 
consider the metric variational approach of our model @ . We note that the possibility of F" (R) being zero for some 
R has not been completely discounted in the literature. Nevertheless, we will impose F" ^ for all R. 

The next constraint that we impose on is that in the limit R oo, we have F{R)/R 1. Writing 

F{R) — R + e(i?), this implies that e/i? as i? ^ oo. This in turn implies that F'{R) ^ 1 as i? —> oo, from which 
we can deduce that e{R) — > const as i? oo. This condition is imposed to ensure that the modified gravity model 
acts suitably like ACDM in the large R regime, where gravity is well constrained by PPN constraints and the CMB. 

With this bound, we have F'{R) ^ 1 as i? ^ oo. However, in ref. [1^ it was found that for a given F{R) theory 
to be stable to perturbations in the large R regime, we must have F" > as i? ^ oo. This means that F'{R) must 
asymptote to unity from below, since F"{R) > implies that it must be an increasing function. Next, we can use the 
fact that we have imposed F"{R) ^ for all R. This means that there can be no turning points in F'{R), implying 
that if F'{R) asymptotes to unity from below, then it must remain less than unity for all R. However, we have 
also imposed F'{R) > 0. From this we conclude that F'{R) must be an increasing function of R, and must satisfy 
0<F'{R)<1. _ 

Another set of constraints on F{R) are discussed in ref. |2^, where the cosmology of F{R) models was considered. 
In order for a theory to be a viable alternative to GR, it must exhibit an acceptable cosmological history, that is it 
must have suitable radiation, matter and late time accelerating eras. It was found that in order for a particular F{R) 
to have a valid matter era, F{R) must satisfy m(r) « +0 and ^ > — 1 at r = —1, where m(r) and r are given by 
m = RF" /F' and r = —RF'/F. In addition, to obtain a late time accelerating era such as we are currently in, the 
function m must satisfy either; (a) m = — r — 1, (\/3 — l)/2 < m < 1 and dm/dr < —1 or (6) < to < 1 at r = —2. 
If these are not satisfied for any R, then the model will not have an acceptable cosmological history. 



IV. A TRIAL F{R) MODEL 

In this section, we give an example F{R) function which satisfies all of the conditions in the previous section. Given 
that F'{R) must be an increasing function, and lie in the range < F'{R) < 1, an obvious choice is 



F'{R) = ^ [1 + tanh(ai? - b)] , (10) 

where a > and b are free parameters. This function has the appropriate behaviour F'{R) —> 1 as R oo and 
F' ^ as R ^ — oo. Integrating this F'{R) gives us 



F{R) = - + —\og[cosh{aR-b)]+A, (11) 
2 Za 

where A is an integration constant. We are looking for F{R) models that exhibit late time accelerating solutions 
in the absence of a cosmological constant. Therefore to specify A, we impose that F{Q) — 0. Hence we find 
A = — (l/2a) Incosh(c) and 



F{R) 



-R 



2a 



log [cosh(ai?) — tanh(6) sinh(ai?)] , 



(12) 



It is this F{R) that will be studied for the remainder of this paper. 

We will begin by considering some of the general properties of our model. First, we look for vacuum solutions to 
the field equations, that is solutions to Q{Ro) = 0. The behaviour of Q{R) is exhibited in Fig. 1(a) We have found 



that Q{R) for our F{R) generically possesses three zeros, for any a and b > 1.2. This means that our F{R) model 
can have three vacuum states, which are Minkowski space (i? = 0), and two de-Sitter vacua, one of which is stable. 

To see these vacuum states more explicitly, we move to the Einstein frame. A plot of the potential V{<j) for our 
F{R) is given in Figs. [2(^[2(b)j We can clearly see that the potential has two minima, one at ^ = (which 
corresponds to the Minkowski vacuum) and the other at some V = Vq > 0. This second minima corresponds to a 



(a) 



(b) 



FIG. 1: (a) is a plot of the function Q(R). This function has three zeros, which correspond to three vacuum states of our model. 
(b) are the functions m (solid), r (dotted) and dm/dr (dashed), plotted parametrically as functions of R. We can clearly see 
that there exists an R such that r = —2 and m < 1, suggesting a late-time accelerating phase. In addition, we have m ~ as 
r ^ — 1, which is the condition required for an acceptable matter era [2^. We have set a = 2, & = 1.5. 
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FIG. 2: (a) The Einstein frame potential. Note the presence of two minima; one at V = 0, and one at some V = Vb > 0. It 
is this minima (shown in (6)) that corresponds to a stable de Sitter solution in the Jordan frame. (&) is V{o) in the small o 
regime. We have set a = 2, b = 1.5. 



stable de Sitter vacuum solution in the Jordan frame, and hence our model has a stable accelerating solution, in the 
absence of a cosmological constant. This is exactly the behaviour we seek. Finally, the third vacuum in the Jordan 
frame corres pon ds to the maxima in the potential 2(a) and is unstable. 

Following [22, we also briefly consider the functions m(r), r and dm/dr for our model, as given in section Hill We 
plot these functions parametrically as functions of R in Fig. 1 (b) We see that there is a region in which m w and 
dm/dr > — 1 for r w — 1, which corresponds to a matter era. As i? ^ oo, we have r — > — 1 and m — > 0, suggesting that 
the condition for a matter epoch is satisfied asymptotically. Further, as R decreases, there is a region where r = —2, 
at which < m < 1. This region corresponds to the late-time de-Sitter phase. Hence our model connects the points 
(—1,0) and (—2, a < 1) in the (r, m) plane, as is required [3l[ for it to have an acceptable cosmic history. 
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V. NEWTONIAN LIMIT 



We now consider the Newtonian limit of our model. In doing so, we will be able to constrain the free parameters a 
and b. The constant a will be fixed by the following two conditions; we require the de-Sitter vacuum solution to be 
i? = i?o ~ 12iJg, and as we will see shortly, we require F"{Rq) <C F'{Rq). Using these conditions, the constant b will 
not be fixed, but we will find that it must be large. 

To see how a and b are determined, we consider the Newtonian limit of modified gravity theories. The full field 
equations for the action ^ are 



F'iR)R^.u - ^F{R)g^, ~ V^V,F'(i?) + g^,uUF'{R) = — T^,. (13) 

We would like to perform a weak-field expansion of these equations around the de-Sitter minima of our model. In 
doing so, we hope to obtain an acceptable Newtonian limit. An analysis has been performed previously in Ref. [28j . 
and expanding the field equations about a symmetric de-Sitter vacuum gives 



F'{Ro)5R^u~\F"{Ro)Rog^u5R-\F{Ro) [sg^, - \g^.u5^-F" {Rq) ,5R+]^g^,U5R^ - ^ (^T^, - ig^.T^ 

(14)' 

where 



'^^M^ = 2 (VpV^Jg,, + V,V"(55,^) + --Ro%. - -^Ragt.u5g - - f^^VJg. (15) 

In order to obtain an acceptable Newtonian limit from these perturbative equations, we require that F'(Rq) ^ 0(1), 
F{Ro) ~ 0{Ro) ^ 1 and F"{Ro) ^ F'{Ro)x^, where x is a macroscopic distance. If this is the case, then the 
perturbative equations ^ will be approximately given by 



F'iRo)6R^. « ^ (^T^. - -g^.Tj . (16) 

Equation is simply the weak field expansion of General Relativity, from which the standard Newtonian limit is 
recovered. The only difference is the coefficient of the first term in (|16p . which is F'{Rq) 0(1) as opposed to unity. 
Note that to obtain the correct limit we must have F"{Rq) <C F'{Rq); if this were not the case then fourth order 
derivative terms such as V ^SR would be present in p6|) . and it has been shown that these terms would lead to 
strong curvature at all lengthscales [lOl . 

For our model (and generically in F{R) theories), we have F{Ro) ^ 0(i?o) ^ F'{Rq) for Rq ^ Hq, which is the 
first condition to obtain an acceptable Newtonian limit. However, the second condition, F"{Ro) ^ F'{Ro)x'^, is not 
generically satisfied, and we must choose our free parameters to obtain the correct behaviour. From (jlOp . we find 
F"{R) to be 



F"{R) ^%ech^{aR-b), (17) 

and hence to ensure that F"{Ro) is very small, we require aRo — 6^1. We can use this condition to fix a in terms 
of Rq and b. To do so, we note that in the vicinity of the de-Sitter minima, R ^ Rq, we can expand F{R) and F'{R) 
for large a = aR — b, 



F(i?)«i?-i^^ll±^-f le-2" + o(e-4«) (18) 
2a 2a 

F'{R) « 1 - + O (19) 

F"{R) w 2ae-2" + Q (e"*") , (20) 

and hence from the vacuum field equations Q{Ra) = 0, we can use our approximations (|18p to give 
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Ro = l2Hl « ^ log (1 + e^^) , (21) 

from which we find a « log (l + e-^*") . This fixes the free parameter a in terms of Rq and the remaining free 
parameter b. Next, we can constrain b from the condition F"{Ro) <C F'{Ro) w 1. Using this, and expanding F"{R) 
in powers of e~^", we find that b must satisfy 



^=5<l±f^>«2fl,. ,22, 

cosh 

which is satisfied for large &, such that 8be~'^'' <C 2i?o- 

To review our conditions; we are working with a model with two free parameters, a and b, and we have imposed 
two constraints on this model, which are that Rq = 12i7Q <C 1 and F"{Ro) ^ 1. These constraints force us to choose 
a and b such that a « j^fp- log (l + e^*") and 8be^^'^ ^ 2i?o- Together, they imply that 6 S> 1, and hence a « We 
note that we have not been forced to specify the constant 6, other than imposing the condition 6 3> 1. This parameter 
controls how closely our model mimics GR; the larger we set b the more suppressed our corrections exp(— 2q;) become. 

Before moving onto the cosmology of our model, we remark that the above weak-field expansion was performed 
around the de-Sitter vacuum R ~ Rq. It has been noted [131 that for local tests of gravity, the expansion should 
instead be around the local energy density of the solar system, which is much larger than Rq. In a subsequent paper 
we will perform a detailed analysis of the local tests of gravity for our model, but for now we simply expand around 
R = Rq. Performing a weak- field expansion around a solar system background will not significantly modify our 
conclusion; we are simply stating that in order to obtain an acceptable weak-field limit for small i?, we must choose 
the constants a and b such that aR — 6 » 1 is satisfied. 



VI. COSMOLOGICAL EVOLUTION 



Having considered the Newtonian limit of our models, we now evaluate the Friedmann equation, and attempt to 
reconstruct the cosmic history of our models. We will find that this model mimics general relativity very closely 
throughout its cosmic history. 

To show this, we first write down the equivalent Friedmann equation for a general F{R) model, obtained from the 
(0,0) component of the field equations 



- S^F'iR) + I + 3Hdt [F'{R}] ^ ^. (23) 

With our highly non-trivial F{R), this equation corresponds to a complex, non linear, third-order differential equation 
in the scale factor a. However, we can use our approximations (jlSp to write (|23p in a much simplified form. We first 
note that the approximations (jl8[) are valid for the entire cosmic history of this model, because they are applicable 
for large a = aR — b. Since R > Rq throughout the cosmic history, and since a is at its smallest value at R = Rq, 
then given that we have constrained a to be large at i? = i?o then it must be large for all R > Rq. 
Using the approximations (|18p in (|23p . we find the following modified Friedmann equation. 



where we have used the condition that b must be large, and hence log(l -I- exp(26)) « 26. The first two terms are 
what we would expect from standard General Relativity with a fine tuned cosmological constant A — fa Rq. The 
remaining terms on the left hand side are corrections to the Friedmann equation, but they are highly suppressed by 
a factor of e~^" <C 1 for all R > Rq. It follows that our model will not significantly differ from standard GR. The 
exponentially suppressed corrections to general relativity that we have found are small in the current epoch, and 
decrease in significance as we track our model back through its cosmic history. Additionally, we note that we still 
have a fine tuning problem, in that we have two free parameters, and we have chosen them specifically such that the 
combination A = ^ satisfies A <C 1. We stress however that the combination b/2a in (|24p is not a true cosmological 
constant; rather our model mimics standard general relativity with a small vacuum energy. 
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The fine-tuning found in our model is a generic feature of modified gravity theories. Since we demand that F"{R) ^ 
for small R to obtain an acceptable weak field limit, it follows that in the vicinity of the vacuum R = Rq, we must 
have F"{R) « 0, and hence F{R) « cR + d, where c and d are constants. This means that for vacuum R ^ Rq, we 
can expand F{R) in the form F(R) « cR + d, where c and d will be some combination of constants in the model. 
This is exactly what we have done for our F{R) model above; we have expanded F{R) as F{R) ~ R — ^ for i? ^ i?o- 
It follows that to obtain an acceptable vacuum solution Rq = 12Hq ^ 1, wc must fine-tune the parameters in the 
model. 

The fact that our model must necessarily mimic GR with suppressed corrections is a consequence of the constraints 
that we have imposed. Specifically, we have used the fact that to obtain an acceptable weak field expansion, we 
require F"{R) ~ for small R. In addition, we have also assumed that in the large R regime F{R)/R 1, which 
suggests that F"{R) ~ as i? — > oo. Since we are considering models which have F"{R) ^ for all i?, then it follows 
that F"{R) must be small throughout the cosmic history. 

We expect that any F{R) function that satisfies the consistency conditions that we have considered will possess a 
Friedmann equation of the form (|24p . that is standard GR plus A with small correction terms. These correction terms 
will not necessarily be exponentially suppressed as they are for the F{R) considered in this paper, but they must 
be suppressed throughout the expansion history of the Universe regardless of the particular form of F{R). However, 
although these models will mimic GR, this does not imply that they are observationally identical to GR if we consider 
all aspects of gravity. For example, it might be possible to observe small corrections to GR by considering solar system 
tests of gravity or the binary pulsar. In a subsequent paper we hope to find distinct observational signatures of our 
model by analyzing local gravity constraints. 

Finally, we observe that the conclusions that we have reached do not apply if we were to relax any of the constraints 
on F{R), for example by allowing F" = for some R in the past. Hence it might be possible to construct an 
acceptable model that does not obey all of the constraints at all times, and hence will not mimic GR throughout the 
entire cosmic history. In addition, it is not strictly necessary for an acceptable modified gravity model to obey the 
consistency conditions for all R, since they can only be probed for R> Rq. We could have an acceptable model which 
has F"{R) = for some R < Rq va the future. However, there is a possibility of such models evolving to a final state 
containing instabilities. 



VII. CONCLUSION 



We have taken the consistency constraints that exist in the literature for modified gravity models, and used them 
to construct a particular modified gravity function F{R). A study of the model has shown that in order to satisfy all 
constraints, it must act like GR with highly suppressed corrections throughout cosmic history. In addition, we find 
that even though we do not have a true cosmological constant in our setup, we must still fine-tune the parameters in 
our F{R) function to obtain an acceptable late time accelerating epoch. Although we have come to this conclusion 
for our specific model, we expect our results to hold true for a large class of F{R) models which satisfy all of the 
constraints that we have considered. 

While we were completing this work we became aware of Refs. [sil [3^ . which come to similar conclusions. In 
Ref. [3l|, it was found that in order to satisfy local gravity constraints, the function m(r) as defined in section Hill 
must satisfy m{r) ^ 10^^® for R ~ Hq. This places a severe constraint upon F"{R), which must be very small in 
the current epoch of the Universe. The model considered in this paper can satisfy the constraint m{r) <C 10^^* by 
making the free parameter b sufficiently large. However, such a stringent constraint is likely to yield a model that is 
virtually indistinguishable to GR. 

In Ref. [3^, a modified gravity function of the form 



F{R)^R-m' 'If J (25) 

C2{^) +1 

was considered, where m, n, Ci and C2 are free parameters. The model considered in this paper shares many similarities 
with (pS)) : the condition F'{R) « 1 was imposed in Ref. [s^ both at i? ~ Hq and as i? — > oo. This implies that the 
function F(R) can be expanded as 



C2 4 \ R J 



(26) 
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for small ^ ^ 1, which is valid throughout the cosmic history of the model. We have arrived at a very similar 
expansion for the F{R) function considered in this paper, (fT8|) . The difference between the two models is that ((26|) 
contains power law corrections to GR, whereas (fT^ contains exponentially suppressed corrections. 
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